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16. Abstract 


The problem of flutter suppression is treated from an energy point of view, whereby 
the energy dissipated per cycle of oscillations is reduced to a quadratic form involving a 
diagonal matrix of energy eigenvalues which are of aerodynamic origin only. A simplified 
binary bending-torsion flutter of a wing strip is investigated by several systems: leading- 
edge (L.E.) control surface, trailing-edge (T.E.) control surface, and combined leading- 
edge — trailing-edge (L.E.-T.E.) control surfaces. In each case the control surfaces are 
allowed to be driven by a linear sensor, a rotational sensor, and a combined linear- 
rotational sensor system. The results show that the flutter of the wing strip can be sup- 
pressed, a single control law being used, over a wide range of reduced frequencies, irre- 
spective of the mass of the system, its stiffness, center -of-gravity location, elastic axis, 
the mode of vibration, and the Mach number (in the subsonic range). Parameters affecting 
the implementation of the control law, such as power requirements, phase lags, and ampli- 
tude gains, are investigated. The extension of the results to three-dimensional flow cases 
is discussed and the use of the flutter suppressor as a gust alleviator is also considered. 
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FLUTTER SUPPRESSION USING ACTIVE CONTROLS BASED ON 
THE CONCEPT OF AERODYNAMIC ENERGY 


By E. Nissim* 

Langley Research Center 

SUMMARY 

The problem of flutter suppression is treated from an energy point of view, whereby 
the energy dissipated by the system per cycle is reduced to a quadratic form involving a 
diagonal matrix of energy eigenvalues which are of aerodynamic origin only. A simplified 
binary bending-torsion flutter of a wing strip is investigated by several systems: leading- 
edge (L.E.) control surface, trailing-edge (T.E.) control surface, and combined leading- 
edge — trailing-edge (L.E. -T.E.) control surfaces. In each case the control surfaces are 
allowed to be driven by a linear sensor, a rotational sensor, and a combined linear- 
rotational sensor system. The results show that the flutter of the wing strip can be sup- 
pressed, a single control law being used, over a wide range of reduced frequencies irre- 
spective of the mass of the system, its stiffness, center-of-gravity location, elastic axis, 
the mode of vibration, and the Mach number (in the subsonic range). The leading-edge — 
trailing-edge control system driven by linear-rotational sensors system is shown to be 
the most efficient system. Estimates are made for the power required to drive the pro- 
posed leading-edge — trailing-edge system and the parameters affecting the implementa- 
tion of the control law (such as phase lags and amplitude gains) are investigated. A 
broad discussion is also included regarding application of the flutter -suppressing system 
to three-dimensional wings, the parameters affecting the reduction of the number of con- 
trol systems required to suppress flutter, and the use of the flutter suppressors as gust 
alleviators. The analytical developments are given in several appendixes and the results 
obtained are illustrated in figures and charts. 

INTRODUCTION 

The recent technological advances made in the field of control systems and the 
increased reliability of control system components offer a new way of treating the prob- 
lem of flutter instability. This approach consists of a rapidly responding control system 

*NRC-NASA Resident Research Associate, now Associate Professor, Technion, 
Israel Institute of Technology, Haifa, Israel. 



which is actuated by the motion of the main surface and which leads to an appropriate 
deflection of the control surface. In this way, use is made of control-surface aerody- 
namic forces to combat flutter instability. Such a system is referred to as an active 
control system. 

Attempts to use active controls were made in recent years mainly in conjunction 
with problems associated with gust alleviation (the objective being to reduce the rigid- 
body response to turbulence) and with mode stabilization (the objective being to reduce 
structural dynamic response to gusts). 

The study, design, and test of gust-alleviation systems have been explored over a 
number of decades and although they have been shown to be effective and technically 
feasible (refs. 1 to 4) they have lagged in practical application. Mode stabilization, to 
reduce structural response, is akin to the problem of flutter suppression. It has been 
seriously considered only during the past decade, and has culminated with the installa- 
tion of active systems on the B-52 aircraft to control the response of the rigid-body mode 
and one elastic mode (first aft body bending) to gust inputs. (See refs. 5 and 6.) No 
attempt will be made to review the extensive literature in the field of gust alleviation 
with mode stabilization (for excellent reviews, see refs. 7 and 8) and mention will only 
be made of some of the contributions and results which have direct bearing on flutter 
suppression. 

Flutter suppression implies, essentially, the control of structural modes only. 
Therefore, mode stabilization methods are of prime importance in flutter suppression. 
Associated problems can be classified within a few categories as follows: 

(a) The determination of the contribution of the different structural modes to a 

given arbitrary deformation of the aircraft 

(b) The determination of a control law which will actuate the control surfaces in a 

certain fashion, proportional to the participation of the structural mode in 
the given deformation, and give rise to a stable control system. An unstable 
control system can arise, for example, when a signal sensed from one struc- 
tural mode gives rise to a control force which actuates a second structural 
mode and which in turn gives rise to a force which drives the first structural 
mode. This condition is known as control system induced instability and is 
caused by the sensor-force couplings. 

(c) The determination that the control forces arising from the determined control 

law are of the type to stabilize the motion sensed 

(d) The determination that the stabilized system is insensitive to changes in flight 

configuration, altitude, Mach number and mode shapes 


2 


Different approaches and methods were proposed and studied for the solution of 
the problems mentioned in (a), (b), and (c). These methods include multisensor systems 
for modal observation, force-sensor modal decoupling networks (very sensitive to 
changes in flight configuration), or mode discriminating sensors to avoid control system 
induced instability, and aerodynamic damping forces for mode stabilization. Linear 
optimal control theory has recently been used (ref. 9) to provide a systematic way of 
treating the problem of control system instability and to increase the control-force sta- 
bilization effectiveness. However, the control law which is derived by this method, and 
which depends on the location of the sensors, cannot often be attained by real systems. 
Furthermore, it requires as many control surfaces as modes considered. Recently, 
Wykes (ref. 10) suggested the identically located accelerometer and force (ILAF) scheme 
to insure control system stability and made use of control "aerodynamic damping" forces 
(that is, control deflection proportional to linear velocity) for modal stabilization. 

These procedures permit some measure of optimization of a specific system at a 
specific flight condition (that is, flight configuration, altitude, and Mach number) and for 
specifically defined mode shapes. However, the problem common to all mode stabiliza- 
tion control systems is a tendency to be very sensitive to changing flight conditions. 

Thus, an optimized control system at one flight condition may either show degradation or 
even give rise to adverse effects at another flight condition. (See ref. 10.) Furthermore, 
even the optimization procedures used for a specific flight condition may often lead to 
control systems which are difficult to realize in practice. Wykes (ref. 10) has used only 
aerodynamic "damping" forces with the ILAF concept and in so doing he imposed a 
severe limitation on the optimization technique. It should also be remembered that 
extensive literature exists (refs. 11 to 13) showing that addition of positive damping can 
turn a stable elastic system into an unstable one. It is therefore concluded that the 
approaches to mode stabilization do not appear to offer an appropriate basis for flutter- 
suppression investigations and that a new approach is desirable. 

Mention should be made here, prior to any elaboration on a new approach, of two 
experimental works carried out with the objective of increasing the flutter speed by the 
use of active control systems. (See refs. 14 and 15.) It is not believed, however, that 
these experiments were methodically planned, and they therefore resulted in very minor 
changes in flutter speed that were often accompanied by poor controllability. It seems 
that extensive analytical work which will lead to deeper understanding of the basic param- 
eters affecting the suppression of flutter by active controls is required before new exper- 
iments are made. 

The approach taken in this paper is based on aerodynamic energy considerations 
and is formulated in a manner insuring lack of sensitivity to changing flight conditions. 
Conceptual studies of some basic control systems are made and their relative merit 
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determined. Figures are furnished partly to illustrate the results and partly to supply 
information for general preliminary design considerations. 

SYMBOLS 


b 


c k 

d 


f /3> f S 




h,h 


hi,h 2 


hs 

i-P 

k 


semichord length 

element i,j, of control law matrix |c] 
scalar multiplier of [gJ matrix 
representative oscillatory amplitude 

generalized forces in /3 and 6 coordinates and required to establish a 
desired control law 

element i,j, of the control law matrix jcij 
bending displacement, positive in down direction 
bending amplitudes defined by sketch (o) 
h at x = x s 

reduced frequency, ^ 


4 , z 8 J 


distance defined in sketch (B2) 


M 


Mach number 


mj = irpb 


4 


s 


n number of degrees of freedom of elastic system 

P energy dissipated by system per cycle 
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Pi 


power required to drive control surfaces 


Qh>Qo!>( 


’max 


t 

V 


generalized aerodynamic force along h, a, (3, and 6, respectively 

maximum stress at a specified location 
reference length or wing semispan 
time 

flight speed 


Wh.Wol 
Wp,W 6 J 


virtual work along coordinates h, a, p, and 6, respectively 


x downstream coordinate measured from midchord point 

x s value of x at which the sensor is located divided by b 

y,z leading-edge displacements of the control surfaces defined in equations (B9) 

and (BIO) 


a, a 


p,6 


X 

H 

p 

U) 


oscillatory angle of attack of wing, or oscillatory amplitude, positive in 
nose-up direction 

value of oi at x = x s 

leading-edge and trailing-edge control surface deflections, respectively, 
positive directions defined in sketch (Bl) 

eigenvalue of ju] 

ith eigenvalue of Juj 

fluid density 

oscillatory frequency 
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Matrices: 


[A^]’[Al] real and ima S inar y P art ) respectively, of aerodynamic matrix A 
[b] inertia matrix 


Bl 

B3 


b 2 

B4I 


submatrices of JbJ defined by equations (C2) 


jcj control law matrix 

[c o ],[Ci],[c 2 ] control law matrices defined by equation (9) 
E structural stiffness matrix 


El 

e 3 


E 2. 

e 4 


structural stiffness submatrices defined by equations (C2) 


{Vy column matrix of forces 

JFjJ force matrix defined by equation (C4) 

[g] control law matrix 

^G 0 J,^Gij,|^G2j control law matrices defined by equation (10) 

[h] complex aerodynamic matrix defined by equation (B16) 


[=1 

[h 3 


,[h 2 ], 

( h 4] 


[n S i],[n s2 ]4 

[N A i]{Na2] 


submatrices of H defined by equations (C2) 


matrices defined in equations (D9) 


[Qr],[Qi] real and imaginary part, respectively, of energy eigenvector modal matrix 
^q^> complex response vector 
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O) 

value of <^qj> at flutter speed 


{>} 

complex amplitudes of response vector 


{'Jr} ■('■‘l) 

real and imaginary part of <£cLq} 


t u ] 

aerodynamic energy matrix defined by equation (A9) 


N 

diagonal matrix of eigenvalues of [u] 


Or) ’00 

real and imaginary parts, respectively, of generalized energy coordinates 

0) 

eigenvector of [u] 


Other notations: 


II 

absolute value 


r i T T 
LJ or 

transposed matrix 


* 

complex conjugate 


LJ 

row matrix 


O 

column matrix 


Subscripts: 



min 

minimum 


max 

maximum 


opt 

optimum 


r 

reference 


Dots over symbols denote derivatives with respect to time. 
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THE ENERGY APPROACH TO FLUTTER SUPPRESSION 


Basic Requirements for a Flutter-Suppression System 

In determining a new approach to the problem of flutter suppression, the following 
goals were set: 

(1) To provide some analytical means which will insure, during the process of 

optimization, minimum sensitivity of control system effectiveness to flight 
conditions 

(2) To provide a simple analytical model which will yield some insight into the 

mechanism of flutter suppression by active controls and still retain the 
essence of the problem 

The simplest possible flutter model is the bending-torsion flutter of a rigid two- 
dimensional strip. The idea of using such a model, although very attractive, seems to 
present some inherent difficulties. For example, what should the inertia matrix or the 
elastic matrix values be ? What values should be assigned to parameters such as mass 
ratio or elastic-axis location and how could these results be correlated, even within some 
approximation such as aerodynamic strip theory, to a three-dimensional wing? These 
points should be carefully considered before discarding such a simple model. The natural 
question which follows is how do the values of inertia or elastic terms, for example, affect 
the flutter instability ? Clearly, these terms affect the natural frequency and the natural 
mode shapes of the oscillating strip and the variation of these parameters over a wide 
range of values will yield, for each frequency of oscillation, a wide range of mode shapes 
which represent, in essence, different combinations of bending and torsional displace- 
ments. Since the frequency is itself a parameter, this statement implies infinite com- 
binations of frequency and mode shapes. It may therefore be possible to accomplish this 
result by a direct variation rather than indirectly (through inertia or elastic terms). 
Indeed, since a wing is a continuous system, one isolated strip on the wing can vibrate at 
many frequencies and many mode shapes and therefore, this variation of frequencies and 
mode-shape combinations is also essential for the extension of the results to a three- 
dimensional wing. Similarly, the sensitivity to changes in flight conditions reduces in 
essence to the effect of changing the mode-shape — frequency combinations. Hence, the 
requirement that the binary system be insensitive to flight conditions implies that the 
increase in stability of the binary system be independent of the mode shapes over a wide 
range of frequencies. There remains now to determine an appropriate approach which 
will allow an easy variation of frequencies and mode shapes and indicate their effect on 
the stability of the binary flutter system. At this stage, an energy approach seems to be 
the natural outgrowth of these requirements. This is true since one can expect that the 
work done per cycle by the system will depend on the aerodynamic forces (which are 
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functions of the reduced frequency for any given Mach number) and mode shapes only, 
whereas the magnitude and sign of the work done per cycle indicate the state of stability 
of the system. 


The Energy Analysis 

In this section the energy concept is developed (the details are presented in appen- 
dix A) without imposing any restrictions on the size of the system and the nature of the 
aerodynamic forces. The results obtained apply to three-dimensional wings as well as 
to the two-dimensional model described. 

Let 



- 0.2 


B + 7 rpb 4 s(A R + iAj) {q} + [e] {q} 


where, at flutter, 


(i) 


(F) =0 

and o> represents the frequency of oscillation; JbJ, the mass matrix; [a r ] and [Aj;J, the 

real and imaginary parts of the aerodynamic matrix, respectively; [e], the stiffness 
matrix; p, the density of the fluid; s, a reference length; b, a reference semichord 
length; and , the response vector. 

As shown in appendix A, the work P done by the system on its surrounding per 
cycle can be written as (eq. (A7)) 



where (from eq. (A2)) 

(„}= {doK"'- |q R + iqjJ e“"‘ (3) 

The sign of P is of importance in the determination of stability, and therefore, it 
will be advantageous to convert equation (2) to a more convenient form. As shown in 
appendix A (eq. (A15)), P can be reduced to the form 

P = \ T 2 Pb 4 »2s(L! R j|VJ{| R } + [SljO] 0i>) (4) 
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or 


P = i + £*) + x 2 ( £ | 2 + 4) + . . . Xn(i 


*nl«ln + « 2 


In 


(5) 


where is a diagonal matrix of the eigenvalues Xj, necessarily real, of the Hermitian 
matrix (eq. (A9)) 



and where the vectors and are defined by the transformation (eq. (A10)) 

{%}=[% *‘%] Or + l «i> <’) 

The matrix + iQjj is a square modal matrix of the principal eigenvectors. Equa- 
tion (5) contains some very interesting features: 

(1) The work P done per cycle has been brought to a principal quadratic form 
(which involves a diagonal matrix) in terms of the system responses. Hence, the sign of 

will be independent of the values of the responses provided it is possible to render all 
the X terms positive. 

(2) The values of X are determined from the aerodynamic matrices only. They 
depend on the assumed modes but are independent of the response of the system. The 
values of X will, however, depend on the reduced frequency of oscillation and Mach 
number. In this way, complete decoupling has been achieved between the aerodynamic 
inputs and the response of the system. 

It can now be seen that a necessary and sufficient condition for the suppression of 
flutter for arbitrary values of responses and is that all the X terms be posi- 
tive. This requirement which stresses the arbitrary nature of the response of the system 
is a severe one in most cases. This statement is true since the values of the £ terms 
of a particular system (in a specified flight configuration) may be such as to render P 
in equation (5) positive despite the existence of some negative values of X. It is essen- 
tial, however, to introduce this requirement to insure the dissipation of energy under any 
conceivable flight condition. The objective of the following work will be to find a relation 
between the main surface oscillatory displacements to the control surface deflections to 
insure that the smallest value of X (that is, X m i n j of the simplified flutter model 
becomes and remains positive and assumes adequately large values over a wide range of 
reduced frequencies and Mach numbers. 
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It is interesting to note here that several investigators have, in the past, used an 
energy approach to treat various aeroelastic problems (refs. 13, 16, and 17) while they 
restricted the energy expressions to real variables only. It now appears that the use of 
complex notation throughout the analysis, as shown in the appendix, is advantageous in 
that it leads to a diagonalized quadratic expression. 

Control System Law for the Simplified Flutter Model 

The two-dimensional rigid strip shown in sketch (a) has two degrees of freedom, 
translation (bending) and rotation (torsion). Therefore, two control surfaces are required 
in order to maintain precise control of the two degrees of freedom of the main surface. 
The main surface is thus allowed to have both leading-edge (L.E.) and trailing-edge (T.E.) 
control surfaces (the choice of a leading-edge — trailing-edge control system is further 


v 


Sketch (a) 

discussed in a subsequent section of this paper). It is possible to see, at this stage, that 
eight parameters (four for each control surface) exist which permit a general deflection 
of the control surfaces with respect to the motion of the main surface. These deflections 
are: 

(1) Deflection of the T.E. control surface actuated by tne in-phase displacement of 

the main surface 

(2) Deflection of the T.E. control surface actuated by the out-of-phase displacement 

of the main surface 

(3) Deflection of the T.E. control surface actuated by the in-phase twist of the main 

surface 

(4) Deflection of the T.E. control surface actuated by the out-of-phase twist of the 

main surface 

The additional four parameters similarly relate to the actuation of the L.E. control 
surface by the in-phase and out-of-phase (90° phase lead) movements of the main surface. 
This relation can be put in matrix form by using the notation illustrated in sketch (b) , that 
is, 
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( 8 ) 



C 11 


'21 



+ l 


G 11 G 12 


G 21 G 22 



Equation (8) represents a control law, where the C and G terms are essentially the 
gearing ratios between the motion of the main surface and the deflection of the control 
surfaces. 


Undisturbed position 



No attempt is made at this stage to consider the practical feasibility of mechanizing the 
control law. Such considerations are left for a later stage of this work, after gaining 
some insight into the relative importance of the different control parameters. 

It is possible to increase the number of parameters by letting the C and G 
terms in equation (8) be functions of the reduced frequency. An expansion of the C and 
G values in terms of reduced frequency will be of the form 


[c] = [C 0 ] + k[cj + k 2 [c 2 ] + . . . 


( 9 ) 


[G] = [ G o] + kjGj] + k 2 [G 2 ] + . . . 


(10) 


It should, however, be noted that the aerodynamic coefficients of both the main surface 
and the control surfaces are quadratic in k"* and the coefficients vary only slowly with 
the reduced frequency k. Therefore, a quadratic representation of C and G should 
be sufficient to represent all the coefficients of the main surface for all practical pur- 
poses. Hence, equation (8) can be looked upon as a special form of equations (9) and (10) 
and can be expected to yield good results at low reduced frequencies. The matrices £c 2 

and j^G 2 J will affect the aerodynamic derivatives of the main surface mainly through the 
control surface aerodynamic inertia terms, and these values are known to be very small 


12 



in wing flutter at the normal range of frequencies. It can therefore be concluded that 
although equation (8) will be used throughout the present work, one might consider the 
use of the following equations if a high range of reduced frequencies needs to be 
considered: 


[c] = [c 0 ] + k[Cj] 

[g] = [Go] + k[Gl] 


In the following discussion, it will be assumed that 


[c] - [Co] 


and = [G 0 ] . 


( 11 ) 

( 12 ) 


Analysis of the Simplified Flutter Model 

The first step in the analysis is to determine the aerodynamic matrices which lead 
to the formation of the aerodynamic energy matrix £uj defined by equation (6). The treat- 
ment of the L.E. control surface presents no special problems since, effectively, it is 
equivalent to a deflection of a large T.E. control surface. A detailed derivation of the 
two-dimensional generalized aerodynamic forces per unit span is presented in appendix B. 
The aerodynamic forces per unit span have the form (eq. (B17)) 


Qh 


c 

h/b 

Qa 

) = 7Tpb 4 w 2 [H]<; 

a 

|Q|3 

p 

Qe 


6 




(13) 


where j^Hj represents the matrix of aerodynamic coefficients, Q represents generalized 
forces, and the subscript refers to the degree of freedom upon which these forces act. 
Substitution of the control law defined by equation (8) into equation (13) yields 


where 



n 


7rpb 4 u> 2 



+ h 2 _c_ 
+ H 4 "c" 




1 - 

hi ; 

H2 

J 

h 3 ; 

h i 


(14) 
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Hence Hj, H 2 , H3, and H4 are 2x2 submatrices of |h]. The generalized forces 
Qp and Q§ need to be known only when evaluating the power required by the control 
system to establish the control law defined by equation (8). The bending-torsion aero- 
dynamic forces can therefore be written as 



where 



iAj] = [Hj + H 2 C + iH 2 G] 


(16) 


Equation (16) shows how [c] and J^gJ affect the value of the matrices [ArJ and j^Ajj . 
Finally, substitution of the values of |ArJ and j^AjJ into equation (6) yields the required 
matrix [u] which, in turn, yields the energy eigenvalues. 


Mention should be made here that the derivation of j^uj for the simplified model 
implies that the control surfaces represent irreversible systems, that is, they can be 
deflected by the control system only and are not affected by other forces. 


As a first step, [uj was calculated herein by using aerodynamic derivatives perti- 
nent to Mach number M = 0, derived from potential theory solutions of Theodorsen, 
Garrick, Kiissner, and Schwartz, and summarized in references 18 to 21. In a later 
stage, Mach number effects were considered in calculating £uj by solving, numerically, 
Possio's integral equation (ref. 22) with the two-dimensional compressible kernel. An 
extensive comparison with the numerical results appearing in references 23 and 24 has 
shown that, in all cases, the accuracy of the numerical results was within 3 percent. 


Optimization Procedure 

The eigenvalue solution of the simplified flutter model yields two values of A, the 
smaller of which is designated by A m i n . Imagine now that A m i n is plotted against k“ * 
for different values of a particular control parameter, for example, c ij. The optimiza- 
tion procedure was initially programed to determine the smallest value of Amin at each 
value of a specific Cij (or Gij) as shown in sketch (c). Comparison was then made 
between points such as 1, 2, and 3 and the optimized value of Cij was defined as that 
value which produced the largest minimum value of A m i n (that is, Cij = rg in 
sketch (c)). The deficiencies inherent in this definition of optimum values made it 
impossible to implement in practice. 
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C ij =r 3 c i j =r 2 C ij- r l 


^min 


1/k 

Sketch (c) 

Sketch (d) illustrates how one may be misled to choose the value of Cjj = r 2 as opti- 
mum, on the basis of this definition, whereas the overall performance of Cij = r^ is 
superior. 



1/k 

Sketch (d) 


The optimization with respect to the control parameters was therefore redefined to 
yield the values of the control parameters which give rise to maximum area under the 
curve of X m i n against k -1 . This definition has possible theoretical deficiencies, such 
as shown in sketch (e). Practice has shown that this definition is workable and that, 
unlike the curve in the sketch, the curves are well behaved. The optimization procedure 


\nin 




1/k 

Sketch (e) 
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consists of the variation of a single control parameter at a time, all the other control 
parameters being kept constant. The value of the control parameter which yields the 
largest area under the curve of A against k" 1 is then assigned to the parameter and 
another control parameter varied in a similar fashion. This whole procedure is repeated 
until convergence of the values of the control parameters is reached. Experience has 
shown that convergence is reached after 3 to 4 parameter sweeps. It is, however, advis- 
able to accompany each optimization with a plot of the optimized A m i n against k"l 
curve, together with sensitivity curves which indicate the effect of the off-design values 
of [c] and [gJ . 


Data and Scope of Optimization 

The binary flutter model was allowed 20 percent chord L.E. and T.E. control sur- 
faces. No variation was allowed in these values and they were kept constant throughout 
the investigation. Similarly, a reference point for sensing the motion of the main sur- 
face was kept constant at 30 percent of the chord (measured from the leading edge). This 
point (x = -0.4b) simulates the location of the sensors on the wing. If it should be desired 
to relate the results to another sensor location (x = Xgb), then the following transforma- 
tion may be used to modify the optimized values of [cl and [g|* 


a= a s 
h = h s - a(x s + 0.4)b I 


or, in matrix form 



-x s - 0.4 
1 



The optimum control law will then assume the form: 



(17) 


(18) 


(19) 


(20) 
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where 


\C + iG} s = \C + i(3 opt 

The flutter model was optimized o 
(0.0128 Ski 19.5). 

The numerical work which follows was carried out in four stages: 

(1) Stage 1, optimization at M = 0: Optimization of X m i n with respect to £c] and 
[g] at M = 0. Comparative studies of optimum conditions were made for T.E. control 
system, L.E. and T.E. controls, and L.E. control only (optimization carried out, in all 
cases, for both single-sensor control law and a double-sensor control law). 

(2) Stage 2, power requirement: Rough assessment of the power required to drive 
the optimized control system and its sensitivity to the control parameters. 

(3) Stage 3, compressibility effects: Mach number sensitivity of Amin> [c] 0 pt 
and Mopt- 

(4) Stage 4, some practical aspects of the optimized control laws: The practical 
achievement of the control law from both a block-diagram point of view and performance 
point of view, together with off-design sensitivities. 

The work was organized to determine the possibilities offered by the analysis and 
to look for any deficiencies which might make the results impractical. Allowance was 
also made for possible interactions between the different stages requiring some com- 
promising results. Where interactions were light, no reoptimization followed since it 
was felt that slight refinements should be left to a later stage. 

OPTIMIZATION AT ZERO MACH NUMBER 
Ranges of Optimization 

The variations of the control parameters were confined, in this stage, within the 
following ranges, to insure small control deflections: 

-0.5 £ C n ^ 0.5 
-1.0 S C 12 ^ 1.0 
-0.5 ^ C 2 1 ^ 0.5 


1 -x s - 0.4| 
0 1 


( 21 ) 


ver a wide range of reduced frequencies 
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-2.5 S C 2 2 = 2.5 
-0.5 2 Gu S 0.5 
-1.0 S G 12 ^ 1.0 
-0.5 ^ G 21 ^ 0.5 
-2.5 = G 22 = 2.5 

As can be seen, the L.E. control surface was allowed a somewhat smaller range of 
variations than the T.E. control surface. Some of these ranges were changed at subse- 
quent stages, on the basis of the results already obtained. Figure 1 shows the variations 
of X m j n and X max with k - * for a wing strip having no control surfaces ^or 

[c] = [g] = o). It can be seen that X m i n is negative throughout the range of k 

(0.0128 = k = 19.5) and that X max is positive throughout the same range of k. Fig- 
ure 1 represents the basic values of X and any improvements due to control deflections 
should be determined by comparison with this figure. 

Results for Trailing-Edge Control Only 

The optimized values of jcj and £g] for the case where only a T.E. control surface 
was allowed are 


fcl - 

0 

0 " 

|Cj - 

-0.35 

-1.9 

r r i _ 

r* 

0 

0 


L G J - 

0.35 

0 . 1 _ 



The graphical descriptions of these [c] and [g] values are illustrated in sketches (f) 
and (g). 


l . 35 h/b 


.35 h/b 



Sketch (f) Sketch (g) 
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Figure 2 shows the variations of X m in with k"* around these optimum values of [c] 
and [g]. Figure 3 shows a similar variation of X max with k _ l around the optimum 
values andj^Gj. 

The following points emerging from these figures are worth noting: 

(1) The value of (x m i n ) 0 p t is only marginally positive (except at high k values) 

and is highly sensitive to off-design values 

(2) The most sensitive control parameters are C22 and G22 

(3) The values of C22 which improve X m i n cause X max to deteriorate 

(4) There is an optimum value of the parameter G21 which, essentially, gives 

rise to aerodynamic damping. (See fig. 2(c).) This result is consistent with 
the observations of Wykes (ref. 10) concerning the effects of damping 

Figures 4 and 5 are basically identical to figures 2 and 3 except for the magnified 
k“* scale. 


By referring back to sketches (f) and (g), it is interesting to note that the main 
effect of the in-phase deflection of the control surface is to counteract any lift buildup; 
that is, the lift increase due to the angle of attack a is opposed by the forces created by 
the deflection of the T.E. control surface. Furthermore, the out-of-phase deflection of 
the T.E. control surface increases the damping in bending, on one hand, while the tor- 
sional damping is reduced on the other. It can therefore be seen that flutter suppression 
is achieved by both reducing the energy input into the system and increasing the dissipa- 
tion of energy. 


Comparison can also be made between this T.E. control system, which requires 
two sensors (to determine a and h/b), and a T.E. control system which makes use of 
one sensor only. The control law of such a single-sensor T.E. control system may 
assume either of the following forms. 



+ 1 


0 

G21 




(22) 


Figure 2 shows the reduction in X m i n caused by either letting C21 and G21 
approach zero. It can thus be seen that the system having the two sensors is superior to 
any single-sensor system. 
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Results for Leading- and Trailing-Edge Control Surfaces 

The optimized values of [c] and jjG] for this case, where the system was allowed to 
have both L.E. and T.E. control surfaces and was activated by two sensors, are 



Illustration of the physical meaning of these £cj and £g] values is shown in 
sketches (h) and (i). Figure 6 shows the variation of A m i n with k“* around these 
optimum values of [c] and £g] whereas figure 7 shows a similar variation of A max . An 
identical representation but with a magnified k“l scale is shown in figure 8. 



Sketch (h) Sketch (i) 


A study of figures 6 to 8 yields the following points which are worth summarizing: 

(1) The optimum value of A m i n is large and positive over the whole range of k~*. 

(2) The off-design sensitivity of the optimized values of £c] and £g] is greatly 
reduced as compared with the T.E. control surface only. 

(3) The most sensitive control parameters which affect both A m i n and A max (in 
practically opposite fashions) are, as before, C 22 and G22- 

(4) The parameter Gj2 mainly affects A m i n whereas the parameter G 21 
mainly affects A max . 

(5) The parameters Cu, Cj2> C 21 , anc * ^11 have only very small effect on 
both A m i n and A max . 
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(6) A variation of G21 can reduce A max to the point where A max = A m i n . Fur- 
ther variation of G21 may cause the interchange between A max and A m i n - Hence, 
there must be a critical value of G21 below which, the variation of G21 will greatly 
affect A m j n . 

(7) The parameters C-q, Cj2, G^, and G12 have optimum values coinciding 
with the limits of their range of variation; thus, an expanded range might improve the 
results. 

The preceding sketch of the activated (optimized) oscillating strip shows again that 
the in-phase control deflections counteract any lift buildup whereas the out-of-phase 
control deflection provides both bending and torsional damping forces. 

Once again, the effectiveness of a single activating sensor was tested. Figures 9 
and 10 show typical results of a variation around such an optimized system. Note that 
the rotational sensor activates all the important control parameters except G21. A 
rotational sensor only reduces effectively to letting Cji = C21 = Gji = G21 = 0 and thus 
causes (through G21) the interchange between A m i n and A max . Clearly, the two- 
sensor system is the more effective sensing system and shows a very large effect on 
''■min- 


Leading-Edge Control Surface Only 

An investigation into the effectiveness of an optimized L.E. control surface activated 
only by two sensors was made. A typical variation of A m j n around the optimum values 
of £cj and £g], shown in figure 11, indicates that a L.E. control surface, on its own, is 
effectively of very little use. 


Subsequent Considerations 

To cut down the amount of work in the following investigation, and on the basis of 
the results obtained so far, it was decided to consider only the two most promising sys- 
tems: the L.E.-T.E. system and the system comprising the T.E. control only. Further- 
more, each of these systems will always be considered as being activated by two sensors. 

Remarks on the Comparison Between the Two-Controls System 
and the Single-Control System 

It is felt, at this stage, that the use of both L.E. and T.E. control systems requires 
an additional justification. This is true since it may be argued that an unfair comparison 
has been made in the simplified flutter model, where an attempt was made to control a 
two-degrees-of-freedom system by a single control surface (T.E. only) and thus justify 
the L.E.-T.E. control system. It may further be argued that the introduction of two T.E. 
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control surfaces, suitably controlled, may produce some improvements similar to those 
of the L.E.-T.E. control system. Sketch (j) shows two such T.E. control surfaces con- 
trolling the bending-torsion degrees of freedom of the strip. 

On the basis of the aerodynamic strip theory, one can divide the strip shown in 
sketch (j) into two strips each having a single T.E. control surface. Since these two 
strips perform an identical oscillation, any optimization technique, including the energy 
approach, will yield 



identical control laws for these two control surfaces. Hence, the two control surfaces 
effectively oscillate as a single T.E. control surface, the performance of which has 
already been investigated. It can therefore be stated, on the basis of this argument and 
the results presented in this work, that the L.E.-T.E. control system is superior to a T.E. 
control system.* The physical reason for this superiority lies in the fact that the T.E. 
control surface has combined beneficial-detrimental effects on the system. On one hand, 
the T.E. control surface is very powerful in decoupling the aerodynamic nonsymmetrical 
cross terms, but on the other hand, it reduces both the aerodynamic damping in torsion 
and the torsional aerodynamic stiffness term (this latter term will be later shown to affect 
the nonactivated strips on an elastic wing). The optimization results reflect the com- 
promise obtained between these conflicting effects. The L.E. control surface, although 
very ineffective in decoupling the nonsymmetrical aerodynamic terms, provides a means 
of adding aerodynamic damping in pitch and positive torsional aerodynamic stiffness 
terms. Hence, these two control surfaces supplement each other and can yield good 
results if properly used. 

*This result appears to be self-explanatory. It has, however, an additional impor- 
tance in that it points out that two T.E. control systems, when located on a wing, will be 
more efficient if the elastic deformations at the controls locations are widely different. 
Two adjoining T.E. controls can therefore be expected to be relatively inefficient at the 
low frequencies (leading to almost identical elastic deformation at the two adjoining sec- 
tions) whereas at high frequencies they may be expected to be more efficient. 
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POWER REQUIREMENT 


Preliminary Considerations 

In evaluating the power required when activating the control surfaces, some further 
definition of the fluttering system has to be made. To avoid analyzing very specific sys- 
tems where force-response amplitude and phase characteristics are of prime importance, 
it was decided to generalize somewhat the problem by attempting to assess the upper 
bound of the power requirement. This assessment can be made by assuming that the 
force-response-phase relationship is such as to require maximum power. Furthermore, 
only the peak power requirement will be evaluated as compared with average power 
required per cycle. The detailed analysis of the power requirement, presented in appen- 
dix C, shows that the upper power bound Pi jm ax can written (eq. (Cll)) 


where 


Pl,max 



C T + iG T 



h = h 0 e iwt 


a = o^e 


,iu)t ! 


(23) 


(24) 


The notation j | means the modulus of the complex number, and jjFjj denotes the sum of 
the 2x2 inertia, elastic, and aerodynamic matrices which appear in the equations of 
motion relating to the control surfaces. (See eq. (C4).) The [FjJ matrix has inertia 

terms (both structural and aerodynamic) proportional to structural elastic constant 
terms, aerodynamic stiffness terms proportional to velocity square and aerodynamic 

damping terms proportioned to cuV. Initial results have shown that the aerodynamic 
stiffness terms (proportional to V^) contribute most to the power requirement. Hence, 
the value of Pi 5 niax will, approximately, be proportional to cuV^ and to the response 
squared. 


Gust power-spectral input considerations indicate that low frequencies, where both 
the input power and the response are large, should be critical from the point of view of 
the power requirement. Similar conclusions can be reached if a uniform cantilevered 
beam is considered to oscillate to produce a specified maximum stress S max at a speci- 
fied location. The inertia loading is proportional to u^d, where d is a representative 
oscillatory amplitude. Hence, 
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Thus, since 


Pi, max a ud 2 V 2 


and using this relation yields 


l,max 


OZ (jO 


V a , 2 ) 


or 


p i, 


S: 


max 


max y2 


fa) 1 - 


Thus, it is again shown that low frequencies will, in general, be critical when the power 
requirement is considered. 


Data for Power Estimate 

The flutter model was assigned the following values to estimate Pi ima x : 

Wing chord, 3.048 meters (10 feet) 

Velocity, 243.84 meters/sec (800 ft/sec) 

Altitude, sea level 

L.E. control, 20 percent chord 

T.E. control, 20 percent chord 

|h 0 /b| = 0.2 

ho = 0.3048 meter = 1 foot 
|o^| = 0.044 radian = 2.5° 
w = 12 radians/second 


r c i 

0.5 

1.0 

i c j - 

-0.05 

1.7 


r G l _ 

-0.5 

1.0 

L G J- 

0.45 

0.2 


Mass ratio, 4. 
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0.04 

0.04 


-0.018666 

0.050666 



0.010666 0 
0 0.010666 


N=[° °" 

L J 0 0 

o o 
o o 

and where [B 3 J, [B 4 J, [E 3 ], and [E 4 J are defined in appendix C. This flutter example was 
extracted from reference 25. The control surfaces were allowed inertia values appropri- 
ate to flat-plate-type controls. 



Effect of Power Requirement 

Figure 12 shows the variation of the power requirement with the variation of the 
control surface parameters around their optimum value. The upper curve represents the 
total upper bound for the power required to actuate a unit span of the wing. The other two 
curves represent the decomposition of this power into its L.E. control and T.E. control 
constituents. The following conclusions can be determined from figure 12: 

(1) The value of Pj max per unit span is very sensitive to the values of Cji, 

C 2 j, and Gji. If it is remembered that these parameters have negligible effect on 
A m in> one can see that the values of 32.3 kW/m (13.2 hp/ft) can be reduced by about 
21.6 kW/m (8.84 hp/ft) if C u = C 21 = G n = 0. 

(2) All the other parameters, except G 2 i, which mainly affect A max contribute 
little to the power requirement. The parameter C 12 is later shown to be of great 
importance and it is worth noting that it has a very small effect on Pi^max- 

It seems therefore, rather unexpectedly, that all the important parameters (C 12 
and those affecting A m i n ) contribute little to the power requirement whereas the parasitic 
terms which do not affect A m i n and G 21 and which affect A max have a large effect 
on the power requirement. 


25 


Reoptimization and Reevaluation of the Power Requirement With 
Unimportant Control Parameters Set to Zero 

On the basis of the forementioned results, it was decided to reoptimize £c] and [gJ 
by constraining Cn = C 12 = C 21 = Gjj = 0 and by extending the range of G} 2 to 
-1.5 S G 12 = 1.5. 

The optimized values obtained were 

0 0 " 

_0 - 1.8 

~0 1.5" 

0.4 0.2 

— i— 

Figure 13 shows the large reduction in power requirement due to these changes. 

The value of 32.3 kW/m (13.2 hp/ft) was thus reduced to 10.7 kW/m (4.36 hp/ft). 

Some additional 25 percent reduction in power can be obtained by using L.E. aerody- 
namic balance obtained by shifting the hinge by 5 percent chord. Figures 14 and 15 show 
the improvement obtained in both A m i n and A max . This improvement followed the 
increased range of Gf2, which, in turn, was increased because of the elimination of Cji 

and Gn and, for the time being, Ci2- Sensitivity tests around those new values of £c] 
and [G] (the figures are not presented here, except for fig. 16) did not show any new 
features which have not already been seen in previous figures. 

It can therefore be concluded that the results obtained are very efficient in terms of 
power requirement and that these upper bounds are small in value and therefore cannot 
present major problems in practical applications. 

COMPRESSIBILITY EFFECTS 

The simplified flutter model was optimized with respect to j^cj andj^Gj under the 
constraints that Cn = Cj2 = ^21 = ^11 = 0 for M = 0.5, M = 0.7, M = 0.8, M = 0.85, 
and M = 0.9. For the case of T.E. control only, no constraints were imposed. In both 
cases negligible changes in £cj and £gJ were obtained for all Mach numbers; hence, it can 
be concluded that the optimized control law remains invariant with Mach number. 

The curves of A m i n at the optimized values of the control parameters show 
changes with Mach numbers. Figure 16 shows such a variation of a T.E. control system 
whereas figure 17 shows similar variation of the L.E. -T.E. control system. It can be 
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seen that Mach number effect is beneficial for the whole k range of the L.E.-T.E. sys- 
tem, whereas the T.E. system shows some improvement, at the high k range, with 
Mach number increase, and a deterioration at the lower range of reduced frequencies. 

It may be noted that the curves of X against k“l for the M = 0 case were com- 
puted by using a numerical aerodynamic solution. It can be seen that the agreement with 
the closed-form solutions shown earlier is generally good, except for the very high val- 
ues of k. (See fig. 17.) At the very high values of k, the closed-form solution for 
M = 0 yields positive X values throughout the whole k range whereas the numerical 
solution yields negative X m i n values for k > 4.6. In the case of the T.E. system, both 
the numerical and analytical solutions yield negative X m i n values for the high ranges of 
k. (See fig. 16.) Another point of interest lies in the effect of compressibility on X m i n 
at high k values. Both figures 16 and 17 show that at high k, all the M ^ 0 curves 
coincide. It can thus be concluded that the compressibility effect is beneficial throughout 
the whole range of k values. 

SOME PRACTICAL ASPECTS OF THE OPTIMIZED CONTROL LAWS 

Preliminary Discussion 

In forming the block diagrams representing how the optimized control laws can be 
achieved, it was assumed that the control surfaces' hinge moments are parameters which 
vary rapidly with flight conditions and which may show substantial deviations from theo- 
retical values. It was therefore decided to form a control system which does not require 
any knowledge of the hinge moments, provided the power requirement is adequate. Fig- 
ures 18 and 19 show the block diagram of such T.E. and L.E.-T.E. control systems, 
respectively. It can be seen that the feedback loop error function is based on the differ- 
ence between the actual control deflection and the desired ueflection as determined through 
the optimized control laws. Two main problems present themselves at this stage: 

(1) How can the frequency co be determined in order to perform the required 
division shown in the block diagrams and how sensitive will the system be to errors in its 
determination ? 

(2) The feedback control loop is, essentially, a second-order system having its 
resonance frequency and damping coefficients. It can therefore be expected that phase 
lags and amplitude changes will be incurred between the desired and actual control deflec- 
tions and that these phase lags and amplitudes will be a function of the frequency of 
oscillation. 

It is therefore necessary to investigate how sensitive the system is to these varia- 
tions and whether practical systems or system components can meet any ensuing require- 
ments. Figure 20 shows the effect of phase lags for the T.E. control system whereas 
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figure 21 shows similar variation of the L.E.-T.E. control system at M = 0.9. How 
sensitive the system is to phase lags at the low range of reduced frequencies can clearly 
be seen. Figure 21, for example, shows that a large positive value of A m i n can be 
changed into a negative value by a mere 20° of phase lag. On the other hand, at high 
reduced frequencies, the phase lag effect becomes much smaller (see fig. 20) and becomes 
negligible at values of k > 1. This result indicates that the feedback control loop needs 
to be designed to have a high resonance frequency, relative to the lowest structural fre- 
quency, when flutter is expected to occur at the low range of reduced frequencies. Lower 
relative feedback loop resonance frequencies are required when flutter occurs at the high 
ranges of reduced frequencies. The effects of the phase lags were investigated at various 
Mach numbers; in every case a similar pattern of variation as in these figures was shown. 

Amplitude Effects 

Figure 22 shows a typical variation of the amplitude effects due to the feedback loop 
for a L.E.-T.E. control system at M = 0.9. It can be seen, once again, that the low 
reduced frequency range shows a very great sensitivity to amplitude gains whereas some 
beneficial effects can be seen at the high range of reduced frequencies for moderate val- 
ues of amplitude gains. These results correlate with the phase-lag results for the reso- 
nance frequency of the feedback loop. Furthermore, it shows that an intermediate value 
of damping is required. Too high a damping value will lead to both a rapid variation of 
phase lags and reduction in amplitude with frequency increase, whereas too light a 
damping, although yielding small variations in phase lags, will cause large amplitude 
variations. 

The T.E. control system behaves in a similar manner as can be seen from figure 23. 

Sensitivity of Optimized System to the Determination of the Frequency 

The determination of the frequency, through the determination of the period of oscil- 
lation, in a manner shown in figure 24, was suggested by Harley Brixey of the Boeing Co., 
Wichita, Kansas. This method yields directly the 1/co term which serves as a multi- 
plier in the controls block diagrams. For a multifrequency signal, it has a tendency to 
yield an average frequency. The problem therefore reduces to the effect on A m in of 
multiplying the whole [gJ matrix by a constant different from 1. Figures 25 and 26 show 
the results of a constant multiplier of on the optimized A m i n curve for the T.E. and 
L.E.-T.E. control systems, respectively, both at M = 0. The T.E. control system shows 
such a low sensitivity that an error of 200 percent is hardly noticeable in figure 25. At 
very high values of multipliers, a deterioration in Amin? which starts at the high k val- 
ues, creeps well into the low k values. Furthermore, the off-design sensitivity, shown 
in figure 29, becomes very much amplified (the multiplier = 18). For multipliers 
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smaller than 1, the values of A m j, n become proportionally smaller but remain positive. 
(See fig. 25.) For the L.E.-T.E. control system, improvement in A m i n can be seen in 
figure 26 for the lower range of k with a deterioration, which becomes relatively large, 
for large values of Ck and high ranges of reduced frequencies. (See fig. 26.) It can 
therefore be seen, in all cases, that inaccuracies in the determination of 1/co of the 
order of 200 percent to 300 percent (that is, Ck = 2 or Ck = 3) have little effect on the 
results. For Ck < 1, a reduction in A m i n is shown in figure 26 but in all cases, A mi n 
remains positive and still assumes relatively large values. Sensitivity tests for the 
L.E.-T.E. control system around the optimum values of [c] and [g] for a large multiplier 
(Ck = 18) does show a decrease in sensitivity (see fig. 29 as a representative example) 
rather than an increase in sensitivity. Similar results pertinent to M = 0.9 are shown 
in figures 27 and 28 for the T.E. and the L.E.-T.E. control systems. These results of 
the 1/co effect on A m i n shed a new light on the damping terms G and the form of 
the control law. This effect is discussed in the following section. At this stage, however, 
it appears that the optimized control law can be achieved in practice by a careful design. 

Some Further Consideration of the Control Law G Parameters 

The effect of the accuracy of the 1/co term, discussed in the previous section, has 
a direct bearing on both the control law and the mechanism of increasing the damping of 
the motion of the fluttering system. These points are discussed herein. 

A different control law . - It has been shown that the division by co gives rise to 
some problems, which although they can be brought to a relatively simple solution, war- 
rant some investigation into the possibility of eliminating the need to determine co. 

Consider the simplified control law 



and modify it to 



where the dots represent differentiation with respect to time. Equation (25) can be easily 
implemented and does not require the determination of co. Since h/b and o- are 
dimensional, equation (25) should be reduced to a nondimensional form. This nondimen- 
sionalizing is done by dividing the [g] terms by a reference frequency co r which lies 
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within the range of the expected frequencies of the system. Thus, equation (25) can be 
written, for sinusoidal oscillation, as 



or, in a different form, as 



Equation (27) is identical to equation (8) except for the factor u>/o> r . Therefore the effect 
of the variation of C^, as discussed in a previous section, is identical to the effect of 
applying the control law defined by equation (26) and moving to a different Ck curve for 
different frequencies. (See, for example, fig. 26.) It should be noted that under this 
definition, oj/w r can assume values larger or smaller than 1. Figures 25 and 29(a) 
which relate to the T.E. control system show that the deterioration of A m i n coupled 
with the high sensitivity to off-design values does not offer any attraction to such a sys- 
tem. On the other hand, the large improvement in X m i n at low k values for the L.E.- 
T.E. control system for — > 1 and the decrease in the off-design sensitivity makes this 

simplified control law a very attractive proposition. It may be possible to allow the val- 
ues of C and G to be dependent on k through the introduction of filters; thus, one 
could choose a small value of co r and yet maintain positive values of A m i n at the very 
high values of w/wr. Some initial work done in this direction yielded some results which 
proved to be very sensitive to Mach number at high k and high w/wr values. Some 
further work using some more sophisticated programing which will allow one to cope with 
the large number of parameters involved in such an optimization is required. 

The damping increase .- A second point emerging from the l/o> analysis indicates 
that although a maximum value of Amin is obtained by varying a single Gij term, a 
monotonically increasing A m i n can be obtained by varying all the Gij terms by a con- 
stant factor (as in the case of the L.E.-T.E. control system) or a shallow decrease in 
A m in ma y result (as in the case of the T.E. control system). This behavior was not 
detected by the optimization procedure since relatively large parametric increments were 
used there. For k values smaller than 1, the G values contribute mainly to the aero- 
dynamic damping. From existing literature (refs. 11, 12, and 13) an increase in damping, 
in a single degree of freedom only, of a fluttering system may show some improvement 
accompanied by subsequent degradation in stability as damping is further increased. It is 
also known that an increase in damping by a constant factor in all the degrees of freedom 
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always leads to an increase in stability. On the basis of these results, it appears that 
for aT.E. control system, a single sensor like Wykes' (ref. 10) linear accelerometer, 
for example, essentially introduces damping in bending, whereas a two-sensor system is 
more efficient in that it enables one to provide some small amounts of damping in pitch. 
The L.E.-T.E. control system allows the control of damping in both degrees of freedom 
and therefore leads to substantial improvements. 

This correlation of results with known effects of linear damping suggests that the 
effect of linear damping on flutter speed may be explained in terms of the X terms and 
that this energy approach may have a wider application than originally was anticipated. 
Some further possible use of the energy approach will be mentioned in a subsequent 
section. 


GENERAL DISCUSSION AND RESULTS 

The Cross Effects Between Active and Nonactive Strips on a Wing 

An isolated strip has so far been considered with the objective of transforming it 
from an energy-absorbing element into an energy-dissipating element. If the intention is 
to span the wing with such activated strips, then the problem is solved since the whole 
wing turns into a dissipative system. If the intention, however, is to put a very small 
number of such strips in some isolated locations along the span, then care must be taken 
to insure that the nonactivated strips are not adversely affected by the activated ones. In 
other words, the "neighboring" nonactive strips should not absorb more energy and offset 
the beneficial effects arising from the dissipative active strips. With this effect in mind, 
possible sources for such adverse interactions are sought. However, any such source 
must be of aerodynamic nature since these forces are the only additional forces that the 
active strips introduce. From these aerodynamic forces, the aerodynamic damping terms 
and the aerodynamic cross-coupling terms are of prime importance in turning the acti- 
vated strip into a dissipative element. Therefore, the possible effects of the direct aero- 
dynamic stiffness terms on the nonactive strips of the wing remain to be investigated. 
These terms do not appear in the matrix ^U J since they are eliminated while it is being 
formed. It should therefore be made clear that any changes made to those direct aerody- 
namic stiffness terms cannot possibly have any effect on the values of X m i n . 

Inspection of the optimized oscillation of both the L.E.-T.E. and the T.E. control 
systems show an upward deflection of the T.E. control with positive values of a (arising 
because of the negative values of C 22 ) as shown in sketches (k) and (l). This upward 
deflection gives rise to a large nose-up pitching moment that leads to a negative torsional 
aerodynamic stiffness term. This negative stiffness term leads to a reduction of the fre- 
quencies of the whole wing as the airspeed is increased and thus affects the neighboring 
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L.E.-T.E. control system 

Sketch (k) 


T.E. control system 

Sketch ( l ) 


strips in an invariably unfavorable manner. In addition, this condition may result in very 
low reduced frequencies which could not have been anticipated on the basis of the natural 
frequencies and velocity of the system; thus, the system is "pushed" outside the k range 
of optimization. 

In the case of the T.E. system, nothing can be done to produce positive aerodynamic 
torsional stiffness terms. The L.E.-T.E. system, however, offers a simple means of 
counteracting the nose-up moment by an appropriate deflection of the L.E. control. It is 
found that a value of C\2 ~ 5.6 eliminates the nose-up moment about the quarter-chord 
point (aerodynamic center). This value of C12 ~ 5.6 insures that the aerodynamic tor- 
sional stiffness of the oscillation strip can never be negative. For wings having elastic 
axes location behind the quarter-chord point, the counteraction of the T.E. control to the 
lift buildup cancels some existing negative aerodynamic torsional stiffness terms. Swept- 
back wings, however, may show a reduction in their positive aerodynamic torsional stiff- 
ness because of this reduction in lift. 

In order to minimize any possible cross effects between the activated strips and the 
nonactivated strips, C\2 should be assigned a value which will offset the strip nose-up 
pitching moment about the aerodynamic center caused by the negative values of C22 
(around k = 0.5). It is therefore hoped that by not allowing the system to assume nega- 
tive aerodynamic stiffness terms and thus maintaining the full structural stiffnesses at 
all speeds, only a small number of strips will be required to eliminate flutter. The 
rigid-body modes form an exception to this rule since, in the absence of any structural 
stiffness terms in these modes, the aerodynamic terms are the sole contributors to the 
stiffness. The special case of rigid-body modes is discussed in some detail in a subse- 
quent section. 

In regard to the aerodynamic bending stiffness, it should be noted that a positive 
value of C21 gives rise to positive bending stiffness and vice versa for negative C21. 
The L.E.-T.E. optimized control law sets C21 = 0 and thus leaves the bending stiffness 
unchanged. The T.E. control law where C21 = -0.35 incorporates some reduction in the 
bending stiffness of the system. It can however be modified, if necessary, since its effect 
on A m in is small. 
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Suppression of Flutter in Three-Dimensional Wings 

The suppression of flutter of a three-dimensional wing can theoretically be treated 
as a direct extension of the two-dimensional flutter suppression case. This treatment 
can be done at various levels of sophistication ranging from the simplified strip theory 
approach which requires some "engineering judgment" to some more elaborate lifting 
surface theory used in conjunction with a modal approach. These approaches are 
described in this section. 

Consider a wing which flutters at a speed Vp and that it is to be increased to let 
Vp lie outside the range of speeds 0 to Vp. The first step involves the inspection of 
the deformation of the wing at Vp with the objective of determining the spanwise loca- 
tion which exhibits the largest product of b^q2. (See eq. (A7).) For constant- chord 
wings, this location may be around the wing tip whereas for tapered wings it will be fur- 
ther inboard. This latter location will form the center line of a L.E.-T.E. strip of length 
Z, governed by the optimized control law. This strip will not only stop absorbing energy 
but will also dissipate a considerable amount of energy. It is further hoped on the basis 
of the precautions taken while the direct aerodynamic stiffness terms were considered 
that the neighboring strips will not extract more energy from the surroundings and thus 
offset the net energy input into the system. The flutter equations are then solved and the 
new flutter speed Vpj is determined. If Vp^ still lies within the range 0 to Vp, a 
second strip is added. This procedure is repeated until the required speed range is 
cleared. It can be seen that this procedure is possible to accomplish since by spanning 
the wing with such strips, the whole wing is turned into a dissipative system. The 
engineering judgment is required to determine the best locations of the strips to keep 
their number to a minimum. 

A second possible approach, still based on aerodynamic strip theory, involves the 
use of the structural mode shapes to compute the energy matrix £uj. Assume that the 
system is allowed to have m mode shapes and therefore m values of A will be 
obtained. Determine the span l of the L.E.-T.E. control surfaces and optimize A m in 
with respect to both the control parameters and j^G^j (in much the same way as in the 
case of the two-dimensional strip) and the spanwise location. Once the optimum control 
parameters and optimum spanwise location are determined, vary the L.E. control param- 
eter Cj2 to offset the negative torsional aerodynamic stiffness effects due to the opti- 
mized value of C22- Determine the important parameters affecting A m i n and discard 
the parasitic terms in the same way as was done for the two-dimensional strip. Reopti- 
mize A m j n with respect to the control parameters constraining the value of C 12 to 
the specified value and the parasitic control terms to zero, the spanwise location 
remaining unchanged. At this stage, the flutter equations are solved again and the new 
value of V F is checked against Vr. If insufficient increase is obtained, an additional 
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strip is added in an identical way. Note that a maximum number of n/2 strips are 
required to control the n modes and that if n is sufficiently large, the optimized con- 
trol law of the strip must be identical to that for the two-dimensional strip results pre- 
sented. This statement is true since if n is sufficiently large, the mode shapes can be 
combined in such a way to produce zero deflections at all the n/2 strips (which span the 
wing) except for the active strip treated. Since only a single strip is oscillating under 
such modal conditions, the optimum control laws must be identical. In this fashion con- 
siderable labor can be saved since the optimized results of the two-dimensional strip 
(with or without the term co/o o v ) can be used right at the start, and the flutter equations 
solved for different spanwise control locations to determine the best location. Alterna- 
tively, the spanwise station which produces the maximum value of Amin is chosen to be 
the best location. 

A third approach is in essence identical to the modal approach which does not make 
use of the preoptimized strip results. The only difference lies in the fact that lifting sur- 
face aerodynamics, incorporating control forces, are used in evaluating the matrix j^uj. 

The strip's values for the control parameters can be used, in this case, to provide the 
range of values around which the optimum is expected to lie. 

It should be stressed here that the number of mode shapes that have to be considered 
in the analysis depends on the design of the control system. The rule should be that all 
the modes that actuate the control surfaces have to be considered. It has already been 
shown (when the c t>/a> r effect was discussed) that a control law which proves to be bene- 
ficial at the low range of reduced frequencies may become detrimental at the high ranges 
of k values. It will, therefore, be incorrect to pursue a flutter analysis involving four 
modes, for example, on the basis that these are the modes which yield flutter within the 
range 0 to Vr. It may well be that the control law may deteriorate the higher frequency 
modes to the extent of producing a flutter speed, at a high k value, within the range 0 
to Vr. It is therefore extremely important to make sure that all the modes which lie 
within the frequency band that actuates the controls are considered. This requirement 
gives rise to a new problem which relates to the practical implementation of the frequency 
bounds of the control law. The lower frequency bound can be established by means of a 
"washout" filter, and no special difficulties arise. The upper bound does give rise to 
some problems since the introduction of a filter for this purpose is accompanied by large 
phase lags and amplitude changes, the effect of which has already been shown to be detri- 
mental. It is therefore felt that a very good way of filtering out the high frequencies is 
through integrations of signals; it should be remembered that the high frequency signals 
are of less importance because of the integrations of the accelerometer signals. If fur- 
ther reduction of the high frequency signals is required, additional integrations can be 
performed on the two signals which form the feedback loop error function. This proce- 
dure is permissible since the control law (eq. (8)) will still be valid if both sides of the 
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equations are integrated an equal number of times. In this way the proper bounds are 
set to embrace the full range of the desired frequencies, and to reduce the high frequency 
signals to the level of the actuator's friction. 

The optimized control law for the L.E.-T.E. control system was applied, analyti- 
cally, on a large delta-wing-type transport airplane at M = 0.9, under NASA contract 
NAS1-9808. The results showed that a single L.E.-T.E. control system covering 12.5 per- 
cent of the wing semispan s and located about 0.1s inboard of the wing tip gave rise to 
28 percent increase in flutter speed. An additional adjoining L.E.-T.E. strip covering 
12.5 percent span and located further inboard led to a complete suppression of flutter in 
the range of speeds that were considered. (That is, the maximum speed tried was 41 per- 
cent greater than the flutter speed of the nonactivated wing. The flutter speed with two 
active strips appears to be much above this speed, if at all existent.) This result is very 
encouraging in that it shows that two L.E.-T.E. strips are sufficient to suppress flutter in 
a real system. 


Suppression of Rigid-Body Modes 

The suppression of rigid-body modes does not warrant a special treatment from the 
point of view of X m i n , provided the reduced frequency lies within the range of optimiza- 
tion and the phugoid motion is not considered. The reservation concerning the phugoid 
originates from the energy approach used in this work. It should be noted that in devel- 
oping the work done per cycle by the system, it has been assumed that the aerodynamic 
forces are linearly dependent on the attitude of the components of the system. This state- 
ment is true for all elastic vibrating modes and it presents a good approximation for all 
the oscillating rigid-body modes with the exception of the phugoid. The phugoid will not 
be treated here and the parameters controlling its motion can be readily identified in the 
control law. With these reservations understood, any strip fitted with the L.E.-T.E. con- 
trol system will dissipate energy when the airplane moves in a rigid-body mode. Special 
consideration, however, may have to be made regarding the contribution of the neighboring 
strips to the energy equation (the assumption being made that the activated strips do not 
span the wing). For illustrative purposes, consider two types of airplanes: one with 
sweptback wings and one with "straight" wings. 

A wing strip with a L.E.-T.E. control system located around the tip of a sweptback 
wing is shown in sketch (m). A positive increase in the angle of attack a, which is nor- 
mally accompanied by an upward increment in lifting force, activates the T.E. control 
upward to counteract the lift buildup. This reduction in lift due to the activated control 
deflection introduces a destabilizing moment about the center of gravity of the airplane, 
or effectively, a relatively large and negative contribution to the rigid-body aerodynamic 
pitching stiffness (due to the large value of Ax). As already mentioned, since the 
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Sketch (m) 


aerodynamic stiffness is the only stiffness in pitch, changes leading to its deterioration 
may affect the response of the rigid airplane in a way which causes the "neighboring" 
strips to absorb more energy and thus, lead to a net destabilizing effect. It can therefore 
be seen that a small reduction in lift at a section located far from the center-of-gravity 
location may cause a large reduction in pitching stiffness. This effect can be readily 
overcome by the introduction of a linear stiffness term due to C21 (which was shown to 
have a negligible effect on the suppression of flutter). It can thus be seen that a small 
bending stiffness term (that is, positive C21) introduced at an activated strip located at a 
large distance behind the center of gravity compensates for this deterioration in stiffness. 
In nonswept wings, this reduction in lift leads to a negligible change in the rigid-body 
aerodynamic pitching stiffness. A review of these considerations which essentially deal 
with avoiding the degradation of the rigid-body modes due to the activation of an isolated 
strip indicates that this problem is of minor importance (since the washout filter required 
to suppress the very low frequency range and allow control of the airplane can be designed 
to filter out the rigid-body frequencies) as compared with the problem which may arise if 
an attempt is made to control the rigid-body modes to improve their stability. 

Assume that it is desired to control, for example, the pitching mode of an elastic 
aircraft. On the basis of the foregoing results, a L.E.-T.E. control strip at the root of 
the horizontal tail (for minimum effect on its neighboring strips) will insure the increase 
in dissipation of the elastic modes of the horizontal tail. Adjustment* of C21 will then 
provide the stiffness control of the rigid-body stiffness in pitch (together with fuselage 
bending stiffness). Any possible adverse effects on the flutter speed due to possible inter- 
action with the wing may be treated by fitting activated strips on the wing at some appro- 
priate spanwise location in a manner identical to that already described. A canard sur- 
face can be treated similar to the horizontal surface. The yaw control also reduces to a 
similar procedure applied to the vertical-tail stabilizer. Finally, it should be noted that 

*When C21 is given a value different from zero, Cji may be introduced. The 
reason for this effect is identical to the one which led to the introduction of Cj 2 due to 
the existence of C 22 - The ratio of C21 should therefore be approximately the 

same as C 12/^22* 
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such active strips, located at the tail surfaces, can be used to compensate any rigid-mode 
degradation due to the wing's active strips, fitted with the objective of suppressing flutter, 
such as the case of the sweptback wing previously mentioned. 

Gust-Response Problems 

The control of the rigid-body modes, as discussed in the preceding section, does 
not give rise to requirements which are contradictory to those necessary to control the 
elastic modes. Hence, no difficulties should arise when attempting to control, simulta- 
neously, both the rigid-body and the elastic modes. Since this type of problem falls 
within the province of gust-response problems, it is evident that the L.E.-T.E. control 
system, properly adjusted for rigid-body control, should also prove to be a very efficient 
gust alleviator. 

As an example, consider the response of an airplane to a vertical gust. It is known 
that the gust problems are usually associated with the rigid-body translation of the air- 
plane, whereas others are associated with its elastic response. A L.E.-T.E. control sys- 
tem located near the root of the wing should lead to the dissipation of large amounts of 
energy due to rigid-body movements. Additional control of rigid-body translation move- 
ment can be achieved, if necessary, by introducing aerodynamic translational stiffness 
through positive values of C21- As noted in the preceding section, the existence of 
C21 ^ 0 may lead to the introduction of Cjj. The problems arising from the elasticity 
of the airplane are best tackled by fitting a L.E.-T.E. control system at a location near 
the tip of the wing n and C21 may again be introduced if necessary^. It should, 
however, be remembered that although Cji and C21 have negligible effects on the 
flutter-suppression characteristics of the system, the penalty they introduce in terms of 
power requirements is by no means negligible. 

Effect of Chordwise Deformations 

The results obtained so far relate to wings with rigid chords. The extension of the 
analysis to a chordwise flexible strip can be made at various degrees of sophistication. 

In the following discussion, only the simplest possible extension is indicated. It is 
assumed that the chordwise deformation can be approximated by three straight lines 
representing the L.E. control surface, the main wing surface, and the T.E. control sur- 
face. This approximation is illustrated in sketch (n). The problem thus reduces to 
determining the control deflection due to the deformation and comparing that deflection, 
in a feedback loop, with the desired deflection determined through the control law. This 
result can be achieved in a simple manner by employing four sensors: two linear accel- 
erometers and two rotational accelerometers. The two linear accelerometers were 
located near the control surface hinges at points 1 and 2 indicated in sketch (o). The 
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Sketch (o) 


output of these two accelerometers can be processed to give a and h/b of the main 
part of the surface. The two rotational accelerometers, are located on the control sur- 
faces, somewhere near the hinge line, at points 3 and 4. (See sketch (o).) Now 




h l ' h 2 
AXJ2 


0e = “4 * “e 
5 e = « 3 - Qfe 

where the subscript e refers to the equivalent deformed chord deflections. (See 
sketch (p).) From the knowledge of and h/b, the signals are processed as before to 
implement the control law, and then the values of the existing control deflections /3 e and 
5 e are subtracted while the error function is formed in the feedback loop. 
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Sketch (p) 

V-g Plots of a Two-Dimensional Flutter System 

The results obtained were applied to a strip flutter problem, the objective being to 
check the control-law effectiveness in terms of a V-g plot. Figure 30(a) shows a V-g plot 
of the binary system previously used for the estimation of the power requirements (taken 
from ref. 25). The additional inertia and elastic data are 

"l.O 0.2 

Bi = 

0.2 0.25 



fo.0625 





the mass ratio is 4 and Bi, B2, Ej, and E 2 are the submatrices used in appendix C. 
The system has a flutter speed Vp/cooto of 1.54 and a divergence speed Vd/wqJd of 2.2 

Figure 30(b) shows the same system activated by the suggested control law; that is, 
C 
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1.5 

0 . 1 _ 

It can be seen that flutter has completely been suppressed and that the g values are in 
the order of hundreds. 

Figure 30(c) shows the same system but with a very small activation that amounts 
to 10 percent of the suggested value. It is interesting to note the resulting major change 
in the V-g plot and the fairly substantial change in flutter speed. This particular result 
indicates that wing designs having large masses, like engines, attached near the wing 
trailing edge may introduce camber changes equivalent to control deflections and may 
lead to flutter results highly sensitive to the mode-shape definition. 

Some Further Use of the Energy Approach 

It has already been mentioned that the energy approach may be applied to investigate 
the effect of linear damping (structural, viscous or aerodynamic) on the flutter stability in 
terms of energy eigenvalues and energy mode shapes. It can also be used in either pre- 
liminary or advanced design to determine those mode shapes which absorb energy from 
the surroundings. 



Determination of the Energy Absorbing Modes 


For preliminary design, it is possible to represent the mode shapes of the structure 
in terms of either synthetic modes or some other known modes taken from a similar air- 
craft. Determine the matrix |u] (eq. (A9)) and solve for the eigenvalues (by using 

eq. (A8)) to form the matrix £aJ and the modal matrices [^Qr + iQij. Determine next, the 
smallest value of X (generally negative) and isolate its corresponding modal columns 
{Tr + iqj}. Inspection of the elements of both this q^ and qj will yield the physical 

modes which participate most in this energy mode shape. If this term is the only nega- 
tive X term, design parameters like engine location, external stores locations, and so 
forth, can be used in a manner which will suppress those modes which contribute most to 
the energy mode. If there are few negative X values this procedure can be repeated and 
preference can be given to those modes which correspond to the lower frequencies 
(smaller number of nodes) and higher negative values of X. 


For an advanced design stage where the flutter mode of vibration {q^} is known, 
equation (A10) can be used to determine the energy vector + i££) at flutter. Substi- 


tution into the energy equation (eq. (A15)) yields the largest negative product of 
^■i(^Rj + ^l\ ^his it* 1 energy mode can then be identified as outlined previously. 
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In a similar manner, the reduction of the larger flutter system into a smaller sys- 
tem can be achieved by ignoring all those modes which either give rise to very small val- 
ues of X (preliminary design) or very small values of Xjfi;^ + (advanced design). 


Mechanism of Flutter 


An additional use of the energy equation is to further the insight into the mechanism 
of flutter. It is shown in appendix D that the flutter equation 


-LO* 


B + 77 pb 4 w 2 s(A R + iAj) + £eJJ ^q} = 0 


can be reduced to the following equation (eq. (DIO)) by the use of the energy modal matrix 
[Qr + iQjJ as a transformation matrix; that is, 

,2 


U 2 [N S1 ] + [n s2 ] + i 


m.w 


— A - cj2n a1 + N A2 


|{«R+i*l} =° 


where the matrices Ngj and Ng 2 are symmetric whereas N A j and N A2 are anti- 
symmetric. All the N matrices and mj are defined in appendix D. The form of equa- 
tion (DIO) is very interesting since it shows that the A matrix (which is diagonal) is the 
equivalent damping matrix. The matrices N A i and N A2 are nondissipative matrices 
since they are antisymmetric and their effect is to provide out-of-phase coupling terms 
only in the equations of motion. Hence, any nonconservative linear problem can be 
reduced in a similar fashion, to a problem of considering on one hand the existence or 
nonexistence of some basically unstable mode shapes (negative X terms) and on the 
other hand, the amount of coupling between those unstable modes and the other stable 
ones. It can therefore be seen that the key point for flutter to be at all possible is the 
existence of a mode shape with negative X. 

The relative importance of the low frequency modes over the high frequency modes, 
when flutter problems are considered, is widely accepted in aeroelasticity (excluding panel 
flutter). The analysis of an elastic system is invariably performed by considering its 
first few elastic modes. The normal justification, based on the responses of the system, 
is not convincing especially when dynamic stability problems are considered. It is, how- 
ever, well known that practice provides the best justification for this approach. Inspec- 
tion of figures 1 and 2 seems to provide some explanation and justification. At a certain 
value of flight speed, the low frequency modes will give rise to low reduced frequencies 
whereas the high frequency modes give rise to high values of reduced frequency. This 
statement means that the low frequency modes will lie on the far right-hand side of fig- 
ures 1 and 2, whereas the high frequency modes will lie on the near left-hand side of 
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these figures. Since A m i n decreases in absolute value continuously and rapidly as k 
increases, it reaches a value of k at which A m i n is completely canceled by the struc- 
tural damping of the system and thus makes it impossible for flutter to exist. 

CONCLUDING REMARKS 

The advantages of using a leading-edge — trailing-edge control system, driven by 
two sensors, and the usefulness of the energy approach have been demonstrated in the 
present work. The results obtained show that flutter can be suppressed but that great 
care must be exercised. It is appropriate to stress here, once again, some of the salient 
points which were brought to light in this work: 

(1) Oscillatory aerodynamic derivatives are essential in this type of work. Fur- 
thermore, they must have good accuracy up to high values of the reduced frequency k. 
The accuracy of the aerodynamic derivatives used in the present work was not sufficient 
at the very high k values, as shown by comparison with computed exact results. Com- 
parisons made with tabulated results, and which showed good agreement, were for values 
of k < 1 (since no higher values of k were tabulated). The results presented in this 
work, using numerical aerodynamic solutions, lead to somewhat smaller values of the 
energy eigenvalue A in the high k range than the exact solution. 

(2) Mach number effects are of a complex nature and no overall correction, like the 
Prandtl-Glauert correction, can be considered satisfactory. An example which showed 
excellent stability from the A m in point of view, over an extremely large range of k, at 
a zero Mach number showed some very adverse effects at a Mach number of 0.9 and 

k > 0.5. Hence, the oscillatory Mach number effects must be considered. 

(3) Control surface and actuator dynamics are of prime importance. Phase lags 
and amplitude changes can easily deprive the system of its effectiveness, especially at 
the low k range. 

(4) All the structural modes that give rise to signals that actuate the control sur- 
faces must be considered in the flutter analysis with active control. This statement is 
true irrespective of the number of structural modes required to represent the important 
flutter modes of the unactivated system. 

It is believed that the control system which consists of the leading-edge — trailing- 
edge controls and which utilizes the frequency ratio co/o> r term (and obviates the need 
to determine the frequencies) is the best from the point of view of practical implementa- 
tion and performance. This condition is particularly true since the very low reduced fre- 
quencies lying around the phugoid or short-period modes do not need such large positive 
values of A m i n as compared with those at the higher values of k. Further 
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improvements in the control law can be obtained if the k range of the system is more 
limited than the one treated here. 

The implications of a successful flutter suppressor on design philosophy are many. 
They span a very wide range of problems: elastic stability, structural efficiency, struc- 
tural fatigue life, and riding qualities of the aircraft. Additional work is, however, 
required before such a system can be incorporated into a design. Analytical work on a 
typical airplane can provide answers to questions relating to the number and size of sur- 
faces required to suppress flutter for different types of airplanes and different categories 
of planforms. As a direct extension of this work, systematic analyses can be done on 
simple types of airplanes where some form of strip theory can be used. The extension 
using lifting surface theories requires the refinement of present numerical schemes to 
insure sufficient accuracy to allow partial span control at high Mach numbers and high k 
values. Finally, some carefully designed experimental work is required to provide sup- 
port to these analytical results. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va. , January 27, 1971. 
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APPENDIX A 


THE ENERGY ANALYSIS 


Let 


{f} = ^-0)2 B + 7rpb 4 s^A R + iAjj + (qj. (Al) 

where, at flutter 

{F> =0 

and co represents the frequency of oscillation; £bJ, the mass matrix; j^A R J and J^Ajj , the 
real and imaginary parts of the aerodynamic matrix, respectively; EJ, the stiffness 
matrix; p, the density of the surrounding fluid; b, the half of a reference chord; {q}, 
the response vector; and s, the semispan of the wing. 

The vector ^q} can be written as 


{q}= {q 0 } e iwt = {q R + iqi}e : 


i oj>t 


(A2) 


where ^q^ is a complex vector of amplitudes. The real part of equation (Al) is essen- 
tially the same as the imaginary part, except for the initial conditions at t = 0. In a 
harmonic motion, where no transients are treated, the real part only of equation (Al) can 
be considered without losing the generality of the equations. Hence, the forces which 
the system exerts on its surroundings are given by 



l-u>‘ 


B + 7rpb 4 s^A R + iA ij + [E]j{q o } 


dcot 




(A3) 


where the asterisk denotes the conjugate vector. The vector {V} is clearly a real 
vector. 

The velocity vector ^q^> can be obtained by differentiating equation (A2) to obtain 


{q} = iw (q 0 ) eiwt 


(A4) 
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APPENDIX A - Continued 


The real part of this velocity vector R^q} is given by 


R {q} = \ {q 0 } eiwt - | iw {q*}e“ iwt (AE 

Hence, the rate at which the system does work on its surroundings is given by R^qJ (f) 
and using equations (A3) and (A5) yields 


R[qJ{F> = tN (" w2 b + 7rpb4g ( AR + iA ijJ + [ E ]j{%) e2: 

" Tl q °j(" w2 b + 7rpb4s ( A R + iA i) + [ E fj {%} 

+ Tl q °j(’ w2 B + 7rpb4s ( A R - iA i)l + [ E ]Wo} 


- xl q °j(~ w2 B + 7rpb4s ( A R " iA lj + [ E ]j ( q o) 


Hence, the work P done per cycle by the system on its surroundings can be found by 
integrating equation (A6) between t = 0 to t=^-. Thus, 


P = l q oj h 2 B + 7rpb4s ( A R " iA l) + [ E ]J{ q o) 


frf- 2 B + ^Pb 4 s( A R + iA^ + [e] ){q 0 ) 


This equation is a scalar equation and therefore the first expression on the right-hand 
side of this equation can be transposed to obtain ^b] and being assumed to be 
symmetric^ 

p . l!£^ q y[.( Al + A T) + i(Ar - {qo > (A7 

Note that the matrices within the square brackets form a Hermitian matrix. 

Now determine the eigenvalues and eigenvectors of the following Hermitian matrix 
[if] extracted from equation (A7): 


K») -MW 
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where 



and represent the vector ^q 0 ^> in terms of the eigenvectors of equation (A8), that is, 


{%} = [Or + iQi]{^R + ill} 


(A10) 


where [Or] and are square matrices whose columns are the real and imaginary 
parts of the eigenvectors of equation (A8), and <0; R + i£j)> are generalized modal coordi- 
nates (defined by eq. (A10)) associated with the aerodynamic energy. 

All the eigenvalue and eigenvector solutions of equation (A8) can be expressed in 
the classical form (ref. 26, p. 6) [q r + iQjjj\J = jujjQp + iQjJ. Postmultiplying this 
equation by + i£j} and premultiplying it by |^q*j or |jj R - i ^lJ pR " iQj' Y ields 


l 4 * • 1{ -J 

Ql - iQf 

[Or + iQljAj^R + i h } - |_£r - i^lj 

Qr - i Ql 

-( A i * A f) 



+ i(A R - , 

$ 

[Qr + iQi]{lR + i£i} 


where A is the eigenvalues matrix, that is, 


(All) 


M 


x 2 



(A12) 



The right-hand side of equation (All) is identical to the right-hand side of equation (A7) 
except for the Ti^pb^co^s/2 factor. Hence, the left-hand sides of these two equations 
can be equated to obtain: 



APPENDIX A - Concluded 


T T 

The matrix Q R - iQj represents the modal matrix of the complex left-hand 
(normalized) eigenvectors and therefore (ref. 26, p. 26) 

Qr “ iQj [Or + iQf] (A 13) 

where is the unit matrix. Hence P reduces to the form 

P = A^S L{r . 14i | a -| 0r + Mi> (A14) 

or 

f < A15 > 

The energy input per cycle into the surroundings has thus been reduced to a quadratic 
form in terms of the generalized energy modal coordinates. 
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APPENDIX B 


DERIVATION OF THE GENERALIZED AERODYNAMIC FORCES 
OF A L.E.-T.E. SYSTEM 


Consider the determination of the aerodynamic forces acting on the system shown 
in sketch (Bl) 


Undisturbed position 



and assume that the forces acting on a somewhat similar system (as used in ref. 16) 
described in sketch (B2) are known. Note the difference between the h, 0 !,j 3,6 coor- 
dinates and the h,Q!,/3,6 coordinates 


Undisturbed position 



Sketch (B2) 

The arrows indicate the direction of either positive displacements, distances, or forces. 
The symbols L and M denote the total lift and pitching moment and L is assumed to 
act through the quarter-chord point. The force acting on the aileron-tab combination 
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APPENDIX B - Continued 


through the aileron L.E. point is denoted by P and the aileron-tab moment is denoted 
by T. The force acting on the tab through its L.E. point is R and the tab moment is 
Q. The forces in sketch (B2), their direction, and points of application are identical with 
those of Smilg and Wasserman (ref. 18) and can be determined directly from expressions 
appearing in the reference for M = 0. These expressions summarize the works of 
Theodorsen, Garrick, Kiissner, and Schwartz (refs. 20 and 21). 

The object of the following analysis is to determine the generalized forces acting on 
the h, a, (3, and 6 coordinates. This object will be accomplished by the application 
of the principle of virtual work. The virtual work in the h coordinate is given by 

W h = b|L (Bl) 

and therefore, the generalized force will be given by 

6W h 

Q h = = bL (B2) 

h 6(h/b) 


Similarly, the virtual work W a in the a coordinate is given by 
W a =Ma- 


Hence, the generalized force Q a is 
SW™ 

Q «=i-r M - u > 

The virtual work in the /3 coordinate is given by 

= -^M - T - P£ 2 )|3 + (L - P)Z 4 J3 

and hence 


Qj3 = 



+ 


Lf 4 + T - M 


(B3) 


(B4) 


(B5) 


Similarly, 


W 6 = Q5 - Rf 3 6 
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and hence, 


Q 0 = -5T = Q - m 3 


(B6) 


These equations can be condensed into the following matrix equation: 

r 


r 


- 





-i 


L 

Qh 


b 

0 

0 

0 

0 

0 


M 

Qa 

Q/3 


-h 

4 

1 

-1 

0 

1 

0 

l 2 ' Z 4 

0 

0 

0 

0 

< 

T 

P 

Q 

Qg 


0 

0 

0 

0 

1 

~h 


R 

L J 







— 1 


^ J 


but the forces are of the form 


r *\ 

L 


b 3 


C- 

h/b 

M 

T 

p 

> = 7 rpw 2 

b 4 

b 4 

b 3 


5 

l 1 

z/bj 

Q 

R 

V / 


b 4 

b 3 


6 

y/b 


(B7) 


(B8) 


where [A^] is a 6 x 6 matrix and z 


and y define the hinge locations and are given by 


z = (Z 2 - 4)3 

(B9) 

y = -4 s 

(BIO) 


Therefore, one can write the following equation based on equations (B9) and (BIO) 


"h/b 


r l 

0 

0 

0 

a 


0 

1 

0 

0 

, j3 


0 

0 

1 

0 

z/b 

> = 

0 

0 

(*2 - 4)/b 

0 

6 


0 

0 

0 

1 

y/b 

V. ^ 


0 

0 

0 

1 

CO 

l "o > 




r 

h/b 

a 




(Bll) 
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APPENDIX B - Continued 


Substituting equations (B8) and (Bll) into equation (B7) yields 


r 

Qh 


~\ 

h/b 

Qq; 

< 

|Q/3 

\ - 7rpb 4 co 2 [A{] ^ 

5 

? 

Qs 

J 


6 

L J 


where 


(B12) 


1 


C A i] = 


-hh 

4/b 

0 


0 

1 

-1 

0 


0 0 0 

0 0 0 

1 (l 2 - 4)/b o 

0 0 1 




1 

0 

0 

0 

0 

0 


0 0 0 

10 0 

0 1 0 

0 (i 12 - 4 )/ b o 

0 0 1 

0 0 -l 3 j b 


(B13) 


Comparison between the coordinates h,a,p,6, and h,a,p,6 
relations: 


yields the following 


h = h - 4 a + 4/3 

a = a - fi 


P = P 


5 = 6 


These equations can be expressed as 


c 

h/b 


1 

-Zi/b 

4 /b 

0 


C 

h/b 

5 


0 

1 

-1 

0 

\ 


/3 


0 

0 

1 

0 


P 

6 


0 

0 

0 

1 


6 



_ 



- 




Substituting equation (B14a) into equation (B12) and rearranging finally yields 


(B 14a) 
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APPENDIX B - Concluded 



where 

-4/ b 4/b 0 

1-10 
0 10 
o (z 2 - 4)/b 0 

0 0 1 

0 0 -4/b 

and p5j'*' is the transpose of [d]. 

If 

[«] = [D] T [A3[D] 

then equation (B14b) can be written in the more concise form: 




(B14b) 


(B15) 


(B16) 


(B17) 
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APPENDIX C 


ESTIMATION OF THE POWER REQUIRED TO ACTIVATE 
THE CONTROL SURFACES 


Consider the equations of motion of a L.E.-T.E. wing strip of unit span which is 
allowed to oscillate in bending and torsion in addition to its control-surface deflections; 
that is, 


r 

0 


-a; 2 (B + 7rpb 4 H] + [e]J {cQ. = / 


(Cl) 




where f^ and fg are external forces applied by the control system on the control sur- 
faces and the matrices are of order 4x4. The first algebraic equations in the matrix 
equations (Cl) have already been used in the energy analysis in conjunction with the con- 
trol law defined by equation (2). The last two equations, when used with the control law, 
yield the forces required for the implementation of the control law. 

Denote 


[B] = 




[H] 


Bl 

; b 2 

b 3 

; b 4 

El 

; e 2 ~ 

e 3 

: e 4 

Hi 

: h 2 

i 

h 3 

i 

; h 4 


(C2) 


where all the submatrices are of order 2x2. Substitute the control law in the last two 
equations of equation (Cl) to obtain 
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-u,2 


B 3 + B 4 C + iB 4 G + 7 rpb^H 3 + H 4 C + il^G^ 


+ E3 + E4C + iE^Gj 


fh/bl 


j f m 


« f, 


(C3) 


By denoting 


[F£J = -0)2 B3 + B4C + iB4G + Trpb^Hg + H4C + il^Gj + [E3 + E4C + iE4Gj 


(C4) 


equation (C3) can be written in the form 



(C5a) 


The matrix 



is complex. 


Its real part will be given by 


R 





(C5b) 


where the asterisks denote complex conjugates. The angular velocities of the control 
surfaces are obtained by differentiating equation ( 8 ); that is, 



(C6) 


The maximum instantaneous power requirement Pi >max is given by 



(C7) 
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APPENDIX C - Concluded 


Equation (C6) indicates that the maximum possible value of R /3,6j will be given by 


r|/ 3,6 i * co 
L^’ J/max 


ho/b ,KI 


C T + iG T 


(C8) 


where the notation | j means the modulus of the complex number written within the lines, 
and 


ft l/bl jho/bl . 


,icot 


(C9) 


a 




Similarly, using equation (C5) yields 


R< 


) f 6 


nJ| h o/ b 
Qbl 


(CIO) 


max 


Substitution of equations (C8) and (CIO) into equation (C7) finally yields 


p l,max 


= co 


h o/ b |’| Q '°l |c T + iG T | 


h o/b|l 

K| 


(Cll) 


Equation (Cll) formed the basis upon which the power estimates were made. 
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APPENDIX D 


TRANSFORMATION OF THE FLUTTER EQUATIONS USING 
THE ENERGY MODE SHAPES MATRIX 


The flutter equation 



B + m^Ap 




where 


m 


j = 77Pb 4 S 


(Dl) 


(D2) 


can be transformed to an interesting special form by the use of the energy modal matrix 
[Qr + iQj] as a transformation matrix, that is, 

( q > = [% + i 9l]0R + i «l} (° 3 ) 


Substitute equation (D3) into equation (Dl) and premultiply the resulting equation by 
Qp - iQ^j to obtain 

(Op - B + mjjAp + iAjj" + [e]Vqr + iQ[]{|p + Hi} = 0 (D4) 


Denote 


A R + A R 


[ A l] ' j|ai + 


4 


' a r 


dl 


4 


At - A- 


?] 


(D5) 


where the first matrix on the right-hand side of equations (D5) is symmetric whereas the 
second is antisymmetric. Substituting equations (D5) into equation (D4) and making use of 
equations (All) and (A13) yields 
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lw^m 


-mor + ‘« i> + fe - iq ? 


- 0,2 


b+ ^Ca r + a5 




+ H)[Q R + iQ I ]{S R+ l{ I } .0 


(D6) 


Denote 


mi 

T S i = B + — (Ar + A- 


T S2 = E 


l R + 


t ai = ^r( A i - a ; T 


(D7) 


where Tgj and Tg 2 are symmetric matrices and is an antisymmetric matrix. 

Substituting equations (D7) into equation (D6) yields 




Kl 


Qr ■ 


-u) 2 ^T 01 + iT A1 ^ + T c 


■SI + 1X A1 + A S2 


Pr + i( ^iT Or + i ^i} - 0 (° 8 ) 


Denote 


[»Sl] “ [^R T S1 Q R + Q^Splj " 
= j^R T S2 Q R + Q I ?T S2 Q l] 

[ n ai] = [QrTsPi - Q? t siQr_ 

^A2j = [Qr T S2 Q I " 9F T ! S2 Q R 


^r t ai^i ■ Qi Tt aPr 


q r t ai q r + ^ Tt aPi 


(D9) 


where j^Ng J and j^Nggj are symmetric and and j^N^J are antisymmetric. Rear- 

ranging equation (D8) and making use of the notations defined in equation (D9) yields 


H 2 [ N Sl] + [ N S2] +i 


mic o* 0 

“V A ' 60 N A1 + N A2 


[ 0r + i ^I } " 0 


(DIO) 


It should be noted that all the real square matrices in equation (DIO) are symmetric 
whereas all the imaginary square matrices are antisymmetric, except for [aJ which is 
diagonal. 
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(a) C n . ( b ) ^12- 

Figure 6.- The effects on A m j n of variation of the L.E.-T.E. control system parameters around 

their optimum values. 
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17k 

Figure 9.- The effects on A m i n of variation of C-q around the optimum 
of the L.E.-T.E. system (Cj2 = c 22 = G 12 = G 22 = °) driven by a linear 
sensor only. 
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Figure 10.- The effects on A m i n of valuation of Cj2 around the optimum 
of the L.E.-T.E. system (C-q = C 21 = G^ = G 2 i = 6) driven by a rota- 
tional sensor only. 






Figure 12.- The effects on the power requirement of variation of L.E.-T.E. 
control system parameters around their optimum values, k = 0.075; 

V = 243.84 m/sec. 
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Figure 12.- Continued. 
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Figure 13.- The effects on the power requirements of variation of the 
L.E.-T.E. control system parameters around their optimum values, 
k = 0.075; V = 243.84 m/sec. 
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Figure 13.- Concluded. 
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Figure 14 .- The effects on A m in of variation of Cji around optimum 
of L.E.-T.E. system ^optimum with Cjj = Cj 2 = C21 = = 0 , 

G21 = 0 . 4 ). 
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Figure 15 .- The effects on A max of variation of Cjj around optimum 
of L.E.-T.E. system (optimum with C11 = C10 = C21 = Gjj = 0 , 
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m 

(a) — = 1; wide k“* range. 



Figure 27.- Variation of A m j n with w/u> r using the opti- 
mum values of C and G from the M = 0.9 results 
(T.E. system). 
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Figure 28.- Concluded. 
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(a) Variation of G 22 ; T.E. system. 
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Figure 29.- The effects on A m i n of variation of the control system parameters 
around their optimum values. Cjj = Cj2 * G-q - ~ - 18; M = 0. 



(c) Variation of G 22 ; L.E.-T.E. system. 
Figure 29.- Concluded. 


110 




V 

tXJQ 


(a) Basic binary system; no active controls. 



(b) Binary system with optimized control 
law; L.E.-T.E. system. 

Figure 30.- V-g plots of a two-dimensional, 
two-degrees-of-freedom (bending- 
torsion) system. 
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(c) Binary system with 0.1 of optimized control law values; 
L.E.-T.E. system. 

Figure 30.- Concluded. 
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